Proximity effect-assisted absorption of spin currents in superconductors 
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The injection of pure spin current into superconductors by the dynamics of a ferromagnetic 
contact is studied theoretically. Taking into account suppression of the order parameter at the 
interfaces (inverse proximity effect) and the energy-dependence of spin-flip scattering, we determine 
the temperature-dependent ferromagnetic resonance linewidth broadening. Our results agree with 
recent experiments in Nb|permalloy bilayers [C. Bell et al, arXiv:cond-mat/0702461]. 



PACS numbers: 74.25.Fy, 74.78. Na, 85.75.-d,72.25.-b 

Cooper pairs in conventional superconductors are spin- 
singlet states and therefore cannot carry a spin cur- 
rent. Some aspects of the resilience of the supercon- 
ducting state against spin-current injection have been 
experimentally demonstrated in hybrid ferromagnct- 
supcrconductor spin valves [1], switches [2], and 7r- 
junctions [3]. In these experiments, the spin current 
flow in the superconducting state can only be inferred 
via charge current measurements. This complicates the 
understanding of the spin current flow in superconduc- 
tors. 

Injection of a pure spin current into a superconduc- 
tor has recently been demonstrated by Bell et al. [4] in 
ferromagnet | superconductor structures under ferromag- 
netic resonance (FMR) conditions, in which the process- 
ing magnetization acts as a "spin pump" [5]. The spin 
angular momentum lost by the ferromagnet can be ob- 
served directly in terms of an increased broadening of the 
FMR spectrum. In this Letter we demonstrate theoreti- 
cally that the spin transport thus measured as a function 
of temperature and device/material parameters offers di- 
rect insight into spin-flip relaxation and the inverse prox- 
imity effect in superconductors. Our theory agrees well 
with the recent experimental results [4], and we provide 
suggestions and predictions for future experiments. 

The theoretical challenge of spin-pumping into su- 
perconductors as compared to normal conductors is 
the strong energy dependence of quasiparticle trans- 
port properties around the superconducting energy gap 
[6]. Also, the energy dependent spin- flip scattering rates 
caused by spin-orbit coupling or magnetic impurities dif- 
fer. Experiments that directly probe spin transport, such 
as Ref. 4, therefore provide unique information about the 
spin-flip scattering mechanism. A complicating factor is 
the inverse proximity effect [7] that suppresses the super- 
conducting order parameter close to a metallic interface 
with ferromagnets like Ni, Co, and Fe. The resulting 
spatial dependence of the superconducting gap requires 
solution of the full transport equations in the entire su- 



perconducting layer. The spin currents measured at such 
interfaces therefore serve as probes of superconducting 
correlations in magnetic heterostructures, and the tem- 
perature dependence of the FMR linewidth near and be- 
low the critical temperature can provide a wealth of in- 
formation about spin-flip processes and superconducting 
proximity physics, with potential implications for differ- 
ent areas of mesoscopic physics. 




FIG. 1: The ferromagnet|superconductor|spin reservoir 
(F|S|R) structure. Precession of magnetization m(t) pumps 
spins into S, which can diffuse and dissipate in R. The F|S in- 
terface spin-mixing conductance for spins polarized transverse 
to the magnetization direction is g x and the S|R interface con- 
ductance is g. The superimposed superconducting gap A(x) 
is suppressed close to the interfaces (inverse proximity effect). 

In the following we develop a theory 
of energy-dependent spin pumping at a 
ferromagnet | superconductor interface and the re- 
sulting spectral spin current flow in the superconductor. 
We consider a diffusive metallic hcterostructure con- 
sisting of a superconducting layer (S) of thickness L 
that is sandwiched by a ferromagnet (F) of thickness 
d and a spin reservoir ("spin sink") (R), see Fig. 1. 
The slowly precessing magnetization m(i) emits a spin 
current that is transversely polarized with respect to 
the instantaneous magnetization direction [5]. The spin 
current that flows through S is immediately dissipated 
upon reaching R. R thus increases the sensitivity of the 
experiments to the spin transport properties of S. R 
represents either a cap of an efficient spin-flip scattering 
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material such as Pt or a large reservoir of a high mobility 
metal [5]. We assume sufficient thermal anchoring so 
that heating from absorbed FMR microwave radiation 
can be disregarded. 

The magnetization dynamics is determined by the gen- 
eralized Landau-Lifshitz-Gilbert equation, 
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Here 7 is the gyromagnetic ratio, H c s is the effective 
magnetic field, Ms is the saturation magnetization, and 
V is the volume of the ferromagnet. The intrinsic dis- 
sipation in the bulk ferromagnet is parametrized by the 
Gilbert damping constant Go- I s is the total spin (i.e. 
angular momentum) current generated by the precessing 
ferromagnet. This loss of angular momentum is equiv- 
alent to an interface contribution to the magnetization 
damping and is observable in terms of the enhanced FMR 
linewidth broadening. Our task is to evaluate the effect 
of superconducting correlations on I s . The results can 
be summarized in terms of an effective resistor model for 
the spin transport. We find an energy-dependent spin 
transport resistance of S in series with the spin-mixing 
resistance r = 1/g of the F|S interface in the normal 
state and the conventional resistance r = 1/g of the S|R 
interface. 

To illustrate the physics we first sketch the results for 
m(i) rotating in the xy-plane and in the absence of spin- 
flip scattering (the derivation for the general situation 
will be outlined subsequently). The magnetization then 
emits a time- independent spin current that is polarized 
along the z-axis [5]. The superconducting condensate 
consists of spin-singlet Cooper pairs. A spin current can 
therefore only be carried in S by excited quasiparticlcs. 
Since the low-energy density of quasiparticle states is sup- 
pressed by superconducting correlations, the spin trans- 
port resistivity is enhanced when S is in the supercon- 
ducting state, resulting in reduced spin injection from 
the ferromagnet. The energy-dependent spin resistance 
is governed by a spectral Ohm's law, 



N(0,E) 



A 



N{L,E) 



(2) 



where /?l = l/(hNoD^) is the effective resistivity of the 
superconductor for spin transport in units of e 2 /h, Nq 
is the density of states at the Fermi level in the normal 
state, Dl(x,E) and N(x,E) are the effective spin diffu- 
sion coefficient and the normalized density of state at po- 
sition x and energy E, respectively [6]. At zero tempera- 
ture, the relevant quasiparticle energy E is determined by 
the FMR frequency which is typically /fmr ~ 10 GHz. 
For BCS superconductors /i/fmr/Ao ~ 0.3 K/T c where 
Ao is the bulk zero-temperature energy gap and T c the 
critical temperature of the superconductor. For small- 
angle precession, the effective "rotation" frequency can 



be introduced as / ~ </>/fmr, where <\> is the angle of 
precession. Thus the relevant energy scale for FMR- 
generated excitations is in practice expected to be much 
smaller than /i/fmRi and the characteristic energy of 
pumped electrons is set by the temperature, see Eq. (3) 
below. At the F|S interface N(x = 0,E) « 1 due to 
the inverse proximity effect (see below). R^ s depends on 
temperature through the local gap A(x,T) which deter- 
mines N(x, E) and pl(x, E). The spin current loss of the 
ferromagnet is consistent with the Gilbert phenomenol- 
ogy in terms of an increased damping parameter G. It is 
determined by the spin angular momentum escape rate 
through S and reads 
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where <?l is the ^-factor, /ib is the Bohr magneton, A is 
the sample cross section area, and /fd is the Fermi-Dirac 
distribution function. 

At temperatures 7 1 « T c , A (a;) as a function of the dis- 
tance from the F|S interface approaches the bulk value on 
the scale of the bulk superconducting coherence length 
£0 = \/fiD/2irkBT c . Since the relevant spin resistivity 
Pl(x,E) and thus R^ s are very large for E < A, £ 
sets the penetration length scale for spin current into the 
superconductor. At low temperatures and L > £0 the 
Gilbert damping (3) will therefore be weakly enhanced. 
On the other hand, at T < T c the gap is suppressed 
throughout S and transport channels at energies E > A 
become accessible. and the Gilbert damping then 

approach the normal state values. 

Spin-flip scattering in S dissipates spin current emitted 
from F, and enhances G by suppressing the back-flow of 
spins into the ferromagnet. The spin-flip length in the 
normal state is given by l s t = y/Dr s i, where D is the 
normal state diffusion coefficient. We take spin-flips into 
account that are caused by magnetic impurities as well as 
spin-orbit coupling at impurities in terms of the spin-flip 
rate l/r s f = l/r m + l/r so [6]. The spin-orbit coupling 
respects the symmetry of singlet Cooper pairs, whereas 
the pair-breaking scattering by magnetic impurities sup- 
presses superconductivity and reduces T c . Below T c , the 
spin-flip rates in S depend on energy. For E < A spin- 
flip rates both due to spin-orbit coupling and magnetic 
impurities are suppressed. For T<T C and L > £oj the 
Gilbert damping will therefore be weakly enhanced. On 
the other hand, for E > A the spin-flip rate due to mag- 
netic impurities is enhanced whereas the spin-flip rate 
due to spin-orbit coupling is similar to that in the normal 
state. We therefore predict a non-monotonic temperature 
dependence of the Gilbert damping close to the critical 
temperature when spin-flip is dominated by magnetic im- 
purities. Experimental data indicate that l s { > £0 for 
typical S. l s { = 48 nm and £0 = 13 nm has been reported 
for Nb [1] (which is used in Ref. 4) whereas l s f = 1.1 (im 
and £ = 124 nm for Al [8, 9]. When L < £ spin-flip in S 
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is therefore inefficient since L < £o < l s { in these materi- 
als. We are then allowed to disregard spin-flip scattering 
[5]. On the other hand, when L ^> l s { the spin current 
never reaches R so that G is governed exclusively by spin- 
flip in S for all temperatures. In the interesting regime 
where i s f sa L, the full theoretical treatment sketched in 
the following has to be invoked in order to compute the 
competing effects that determine G. 

The total spin current leaving the ferromagnet in the 
F|S|R heterostructure can be expressed as an energy inte- 
gral over the balance of the spectral pumping and back- 
flow currents I s = / dE(i^ — ig ack ). The spin current 
injected into S by the precessing magnetization is [5, 10]: 

.h™ fi^(0, E) f F p(E - hf/2) - fME + hf/2) 

h {E) hf 

/ I dm I dm\ 
xj^mx— j, (4) 

where / is the instantaneous rotation frequency. Here, 
<7r and g^- are the real and imaginary parts of spin- 
mixing conductance. For metallic interfaces, ^> g^ 
[11]. We therefore disregard the "effective field" in 
(4), although it contributes to the interface boundary 
conditions discussed below. The magnetization damp- 
ing that follows from (4) is frequency dependent beyond 
the Gilbert phenomenology. We have checked numeri- 
cally that the /-dependent terms contribute weakly to 
the damping even when hf < Ao for the parameters 
studied. We therefore restrict attention to the linear 
response regime in which the Fcrmi-Dirac functions in 
(4) can be expanded to first order in hf. This leads to 
frequency- independent enhanced Gilbert damping in (1). 
The spectral back-flow of spin current into F induced by 
the spin accumulation on the S side is 

i h^ {E)= _N(^p_ g ± hTs(0iE) _ (5) 

The nonequilibrium spin distribution function hxs(x,i?) 
can be computed by Keldysh transport theory [6]. 
In the S bulk, the total spin current l s (x) = 
hAN J^° oo dED lj (E,x)d x h.Ts(x, E)/2 follows from the 
diffusion equation 

(Nd t + d x D h d x ^sts _ ^sts \ hxs = (6) 

Diffusion through S is taken to be instantaneous on 
the scale of the FMR frequency as long as / < D/L 2 
and/or / -c l/r s f so that hxs in (6) becomes time- 
independent, a^g^g = [Re cosh 9} 2 + (— )[Re sinh^] 2 are 
energy-dependent renormalization factors for the spin- 
flip rates due to magnetic impurities (spin-orbit cou- 
pling), and the energy dependent spin diffusion coeffi- 
cient Dl/D = a^grpg. The spectral properties of the 
superconductor parametrized by 9(x, E) are determined 



by the Usadel equation for the retarded Green function 
G R = T3 cosh 9 + 1T2 sinh 6, 

h D f) 2 f) 3 h 

— ^— = iAcosh(0) - iE whiff) + -— sinh(20), (7) 
2 ox 2 8 r m 

to be solved with the BCS gap equation A = 
(N a X/2) J a E ° dEtanh(E/2k B T)Rc sinh(0) [6]. Here, E D 
is the Debye cut-off energy and A the interaction param- 
eter. 

The boundary condition for the diffusion equation (6) 
is conservation of spin current at the interfaces. At x = 0, 
hAN D L d x h TS /2 = ifj - i^ ack . We use boundary con- 
ditions derived in Ref. 12 for (7) at the S|R interface. 
At the F|S interface we impose complete suppression of 
superconducting correlations, 6{x = 0, E) = for the 
following reasons. The large exchange energy in transi- 
tion metal ferromagncts completely suppress supercon- 
ducting correlations, so that the F adjacent to S is a 
source of incoherent particles. Additionally, spin de- 
pendent interface scattering at the S side [13] induces 
an effective pair-breaking exchange field, which we es- 
timate as B c tf = hg^-/e 2 gLHBNoA!;o [14]. Here, N A^ 
is the number of states at the Fermi energy within £0 
from the interface. With g^ ss 0.05#sh, where gsh is 
the Sharvin conductance [11], and approximating TVo 
by the free-electron value, //B^eff is comparable to Ao, 
c.g, ^ B B off (Nb) - 0.56 meV, n B B eS (Al) - 69/icV. The 
bulk F exchange splitting and the induced B e s by spin- 
dependent interface scattering leads to a vanishing gap 
(and 9) at the F|S interface [15, 16]. 

The spin diffusion equation (6) can be solved analyti- 
cally in the absence of spin-flip, proving (2). We now use 
the full machinery sketched above to make contact with 
experimental results for a F|S device (without R) similar 
to sample C in Ref. 4. Numerically computing I s includ- 
ing spin-flip caused by magnetic impurities [17], we ob- 
tain the enhanced Gilbert damping G from (1). In the ex- 
periment, F is a permalloy layer with d — 2 nm, and <7l = 
2.1. S is Nb with L = 70 nm, bulk critical temperature 
T c0 = 8.91 K, Z sf = 48 nm, and D = 5.41 cnrV 1 [1, 18]. 
For the interface conductances we use Ar = 3 film 2 [19]. 
We find G - G = 0.777 x 10 8 s" 1 at T c /2 = 3.6 K and 
1.19 x 10 s s _1 in the normal state. When the inhomo- 
geneous linewidth broadening is small, the width of the 
FMR spectra are proportional to G and the experimen- 
tal data gives [G(T > T c ) - G(T = T c /2)]/G(T > T c ) w 
21 %. Using Go = 0.7 x 10 8 s" 1 [5] we obtain 22 %. The 
measured reduction of the Gilbert damping upon cooling 
the sample from above T c to T c /2 agrees quantitatively 
with our calculation. 

We can make additional predictions for the Gilbert 
damping in F|S|R systems, focusing on Al as S since its 
spin-flip length is much larger than that of Nb, and as 
a weak coupling superconductor is better described by 
BCS theory. The Al material parameters are T c o = 1.26 
K, l s { = 1.1 /im, and D = 160 cm 2 s _1 . In the left panel 
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of Fig. 2 we show the temperature dependence of G — Go 
for three different thicknesses L when spin-flip is induced 
exclusively by cither magnetic disorder or spin-orbit cou- 
pling to impurities. In contrast to spin-orbit scatterers, 
magnetic impurities reduce T c due to the pair-breaking 
term in (7). For L > l s { and T <C T c , as well as for 
T > T c , the results do not depend on the nature of 
the spin-flip scattering. In general, we observe that T c 
strongly depends on L due to the inverse proximity ef- 
fect. We also note that the difference in damping between 
the normal state and the superconducting state is small 
when L ~ £ since only a small gap develops. 

The experiments of Ref. 4 probed the regimes L <C £o 
as well as L 3> £o- We also present results for arbitrary 
L/£o- In the normal state, G decreases with increasing 
L due to increasing bulk spin transport resistance, which 
limits relaxation in R, until L reaches the value of l s { 
where R becomes irrelevant (inset Fig. 2). When T<T C , 
on the other hand, the relevant length scale for spin pen- 
etration into S is £o- This explains the more rapid decay 
of G — Go as a function of L in the superconducting state. 
When L > £ , the spin-current absorption is completely 
determined by the inverse proximity effect: Spin dissipa- 
tion in R by transport through S is suppressed by the 
superconducting gap, and, furthermore, spin relaxation 
deep in S is suppressed by the superconductivity. How- 
ever, the inverse proximity effect enhances the density of 
states at low energy as well as spin-flip scattering rates 
close to the F|S interface. 

When L < l s {, the results depend strongly on the S|R 
contact described by g. In the right panel of Fig. 2, we 
show the temperature dependence of G — Go for L = 
900 nm in an F|S system (no R or g = 0). At T > T c , 
the damping is much smaller in the F|S system (the right 
panel) than in the F|S|R system with the same L (the 
middle pair of curves in the left panel). T c is also higher 
since there is no inverse proximity effect at x = L. At 
very low temperatures, T<T C , G — Go saturates at the 
same value for the F|S system as the F|S|R system with 
the larger thickness, L = 1300 nm. For such thick S, 
T c is unaffected by R and spins cannot diffuse through S 
and dissipate in R, so that the resulting damping is the 
same as in the F|S system. We also see from the right 
panel of Fig. 2 that when T < T c the enhanced Gilbert 
damping can be somewhat larger than above T c when 
spin- flip is induced by magnetic impurities, because the 
induced spin accumulation of quasiparticles with energy 
k^T > A experiences an enhanced spin-flip rate through 
a TSTS- I n the F|S|R system, this effect is overwhelmed 
by the spin accumulation drain in R. 

In conclusion, our theory quantitatively repro- 
duces the measured FMR linewidth broadening in 
ferromagnet | superconductor structures. We make addi- 
tional predictions for varying system sizes and temper- 
atures, and the nature and strength of spin-flip scatter- 
ing. We hope to stimulate more experiments that should 
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FIG. 2: Calculated G - G [10 8 s _1 ] (same ordinate in all 
plots). Red solid (green dashed) lines for system where r s f = 
Tm (Tsf = Tso). Left panel: F|S|R system with L [nm] from 
top to bottom: 600, 900, 1300. Right panel: F|S system (no 
R) with L = 900 nm. Inset: L dependence [/tm] of G — Go 
for T > T c (green dashed line) and T <§; T c (red solid line). 



reveal information about the strong inverse proximity ef- 
fect and energy dependence of spin flip scattering in these 
systems. 
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